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ABSTRACT 

Mayall II = Gl is one of the brightest globular clusters belonging to M31, the Andromeda 
galaxy. Our observations with the Wide Field and Planetary camera WFPC2 onboard the 
Hubble Space Telescope (HST) provide photometric data for the I vs. V — I and V vs. V — I 
color-magnitude diagrams. They reach stars with magnitudes fainter than V = 27 mag, with a 
well populated red horizontal branch at about V = 25.3 mag. 

From model fitting, we determine a rather high mean metallicity of [Fe/H] = -0.95 ± 0.09, 
somewhat similar to 47Tucanae. In order to determine our true measurement errors, we have 
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carried out artificial star experiments. We find a larger spread in V — I than can be explained by 
the measurement errors, and we attribute this to an intrinsic metallicity dispersion amongst the 
stars of Gl; this may be the consequence of self-enrichment during the early stellar/dynamical 
evolutionary phases of this cluster. So far, only u Centauri, the giant Galactic globular cluster, 
has been known to exhibit such an intrinsic metallicity dispersion, a phenomenon certainly 
related to the the deep potential wells of these two star clusters. 

We determine, from the same HST/WFPC2 data, the structural parameters of Gl. Its 
surface brightness profile provides its core radius r c = 0.14" = 0.52 pc, its tidal radius r t ~ 54" 
= 200 pc, and its concentration c = log (rt/r c ) ~ 2.5. Such a high concentration indicates the 
probable collapse of the core of Gl. KECK/HIRES observations provide the central velocity 
dispersion <r b s = 25.1 kms -1 , with er p (0) = 27.8 kms -1 once aperture corrected. 

Three estimates of the total mass of this globular cluster can be obtained. The King-model 
mass is Mr = 15 x 10 6 .M© with M/L v ~ 7.5, and the Virial mass is Mwr = 7.3 x 1O 6 7V4 
with M/Ly — 3.6. By using a King-Michie model fitted simultaneously to the surface brightness 
profile and the central velocity dispersion value, mass estimates range from Mkm = 14 x 
10 6 X Q to 17 x 1O 6 A4 . 

Although uncertain, all of these mass estimates make Gl more than twice as massive as 
oj Centauri, the most massive Galactic globular cluster, whose mass is also uncertain by about 
a factor of 2. Gl is not unique in M31: at least 3 other bright globular clusters of this galaxy 
have velocity dispersions a b s larger than 20 kms -1 , implying probably similar large masses. 

Such large masses relate to the metallicity spread whose origin is still unknown (cither 
self-enrichment, an inhomogeneous proto-cluster cloud, or remaining core of a dwarf galaxy). 
Let us consider for Gl (see Table 1) the four following parameters: central surface brightness 
/i(0,V) = 13.47 mag arcsec -2 , core radius r c — 0.52 pc, integrated absolute visual magnitude 
My — -10.94 mag, and central velocity dispersion cr(0) = 28 krns" 1 . When considering the 
positions of Gl in the different diagrams defined by Kormendy (1985) using the above four 
parameters, Gl always appears on the sequence defined by globular clusters, and definitely away 
from the other sequences defined by elliptical galaxies, bulges, and dwarf spheroidal galaxies. 
The same is true for to Centauri. 

Little is known about the positions, in these diagrams, of the nuclei of nucleated dwarf 
elliptical galaxies, which could be the progenitors of the most massive globular clusters. The 
above four parameters are known only for the nucleus of one dwarf elliptical, viz., NGC 205, 
and put this object, in the Kormendy's diagram, close to Gl, right on the sequence of globular 
clusters. This does not prove that all (massive) globular clusters are the remnant cores of 
nucleated dwarf galaxies. 

At the moment, only the anti-correlation of metallicity with age recently observed in 
co Centauri suggests that this cluster enriched itself over a timescale of about 3 Gyr. This 
contradicts the general idea that all the stars in a globular cluster are coeval, and may favor the 
origin of to Centauri as being the remaining core of a larger entity, e.g., of a former nucleated 
dwarf elliptical galaxy. In any case, the very massive globular clusters, by the mere fact that 
their large masses imply complicated stellar and dynamical evolution, may blur the former clear 
(or simplistic) difference between globular clusters and dwarf galaxies. 

Subject headings: galaxies: dwarf, evolution, formation, Local Group, star clusters; Galaxy: 
globular clusters: general Galaxy: globular clusters: individual (u Centauri, Mayall II = Gl 
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1. Introduction 

From both stellar population and stellar dynamics points of view, globular clusters represent a very 
interesting family of stellar systems. They are ancient building blocks of galaxies and contain unique 
information about old stellar populations and their formation. Some fundamental dynamical processes take 
place in their cores on time scales shorter than the age of the universe, offering us unique laboratories for 
learning about two-body relaxation, mass segregation from equipartition of energy, stellar collisions, stellar 
mergers, and core collapse. See Meylan & Heggie (1997) for a general review. 

In our Galaxy, globular clusters span a wide range of properties (Djorgovski & Meylan 1994). For 
example, their integrated absolute magnitudes and total masses range from = -10.1 and hAtot — 5 x 
1O 6 A^ (Meylan etal. 1994, 1995) for the giant Galactic globular cluster wCentauri down to M£?* = -1.7 
and Mtot — 10 3 X Q for the Lilliputian Galactic globular cluster AM-4 (Inman & Carney 1987). AM-4, 
located at ~ 26 kpc from the Galactic centre and at ~ 17 kpc above the Galactic plane, does not belong 
to the Galactic disk, and consequently cannot be considered to be an old open cluster. The uncertainties 
on the above total mass estimates, typically as large as 100%, do not alter the fact that, in our Galaxy, the 
individual masses of globular clusters range over three orders of magnitude. It is not known to what extent 
these mass differences are "congenital" or due to subsequent pruning by dynamical evolution. 

With its approximately 450 members (Barmby et al. 2000), the globular cluster system of M31, 
the Andromeda galaxy, is about three times as rich as the Galactic one, and is among the most studied 
globular cluster systems in external galaxies (Harris 1991). However, our knowledge comes mainly from the 
integrated photometric and/or spectroscopic properties of these clusters (e.g., Reed etal. 1994, Barmby et 
al. 2000). It is essentially since the advent of the Hubble Space Telescope (HST), with its post-refurbishment 
camera WFPC2, that we resolve clearly some of the M31 clusters into individual stars, giving access to 
their general morphology and structural parameters (e.g., Fusi Pecci etal. 1994, Grillmair etal. 1996), 
and providing their Color-Magnitude Diagrams (CMDs), reaching magnitudes fainter than the Horizontal 
Branch (e.g., Ajhar etal. 1996, Fusi Pecci etal. 1996, Rich etal. 1996, Jablonka etal. 2000). Correlations 
between structural, photometric and dynamical parameters have been investigated for 21 globular clusters 
in M31 (Djorgovski etal. 1997). 

The brightest globulars in M31 are brighter than our Galactic champion wCentauri. Among these 
giants is Mayall II = Gl (Rich etal. 1996), a globular cluster so bright that, like wCentauri, it has been 
considered as the possible remaining core of a former dwarf elliptical galaxy which would have lost most 
of it envelope through tidal interaction with its host galaxy (Meylan etal. 1997, 2000, Meylan 2000). We 
present in this paper a detailed photometric and dynamical study of Gl. 

This paper is structured as follows: Section 2 gives some general information about Mayall II = Gl, 
Section 3 describes the observations, Section 4 gives the CMD of Gl and discusses the spread in metallicity 
observed among the stars of the red giant branch, Section 5 presents its ellipticity, position angle, and 
surface brightness profile, Sections 6 and 7 give the results from the various mass estimators. All results are 
summarized and discussed in Section 8. 



2. Mayall II = Gl, a luminous globular cluster in M31 

The globular cluster Gl belongs to our companion galaxy, Andromeda = M31. Resolved with difficulty 
from the ground (Djorgovski 1988, Heasley etal. 1988, Bcndinelli etal. 1990), this huge swarm of stars 
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appears as a bright flattened star cluster when observed with the Wide Field and Planetary camera WFPC2 
onboard the Hubble Space Telescope (HST) (see Fig. 1). The integrated visual magnitude of this cluster, V 
= 13.48 mag, corresponds to an absolute visual magnitude My = -10.94 mag, with E{B — V) = 0.06 and a 
distance modulus (m — M) M31 = 24.43 mag, implying a total luminosity of about L v ~ 2 x 10 6 i Q (Rich 
ctal. 1996, Djorgovski etal. 1997). 

The coordinates of Gl, when compared to the coordinates of the center of M31 (see Table 1), place 
it at a projected distance of about 3°, i.e. 40 kpc, from the center of M31. This rather large projected 
distance is comparable to the distance between our Galaxy and the Large Magellanic Cloud. Nevertheless, 
both color-magnitude diagrams and radial velocities of Gl {V r (Gl) = —331 ± 24 krns" 1 ) and M31 
(K(M31) = -300 ± 4 kms -1 from 21cm HI line and V r (M31) = -295 ± 7 kms -1 from optical lines), 
completely support the idea that this cluster belongs to the globular cluster system of M31. The values of 
the most important general parameters describing Gl are displayed in Table 1. 



3. Observations 

We use our observations obtained with the Planetary Camera (PC) of the HST/WFPC2, in the 
framework of a programme (PI Pascale Jablonka, ID = 5907) aiming at studying star clusters and stellar 
populations in M31 (see Jablonka etal. 1999, Jablonka etal. 2000). The PC pixel scale is 0.045"pix _1 . The 
four images of Gl, taken with each of the F555W (V) and F814W (I) filters, have total integration times 
equal to 500 + 500 + 600 + 600 = 2,200 seconds in V and to 400 + 400 + 500 + 500 = 1,800 seconds in 
/. Because of possible non-linearity in the bright concentrated core of Gl, our rather deep exposures are 
supplemented with some shorter exposures from another programme (PI R. Michael Rich, ID = 5464). See 
Rich etal. (1996). 

Fig. 1 displays an area of 31.5" x 31.5" from the original PC frames centered on the cluster. This 
image is a composite of all our V and / frames and provides a genuine indication of the relative colors of 
the stars. Although completely resolved, the cluster appears extremely compact, with a very steep surface 
brightness profile and an extremely bright and crowded core. 

We determine the photometry using one of the presently best available algorithms for performing 
stellar photometry in crowded fields, viz., the ALLFRAME procedure developed by P. Stetson (1994). 
ALLFRAME is run on 700 x 700-pixel (31.5" x 31.5") sub-areas of the original (36" x 36") PC frames, 
avoiding the egdes of the frame and masking the two areas disturbed by the two bright foreground stars 
(Fig. 1). As ALLFRAME is now widely known and since our use of it is already described in Jablonka et al. 
(1999), here we mainly focus on the results. Given the very large projected distance between Gl and the 
core of M31 (40 kpc), it is worth mentioning that the number of M31 field stars in our PC field is negligible 
when compared to the number of Gl stars. 



4. The Color-Magnitude Diagram of Mayall II = Gl 

Fig. 2 displays the two color-magnitude diagrams (CMDs) of Gl, each of them containing the same 
4903 stars. The left panel shows the V vs. V - I CMD, while the right panel displays the I vs. V - 1 CMD. 
The brightest stars at V~22.5 have color errors of ^0.03, and stars at the level of the HB have color errors 
of ~0.15 mag. 
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A first CMD of Gl, reaching stars below the horizontal branch (HB), was published by Rich etal. 
(1996) based on HST Cycle 4 data, with total exposure times of 1,600 seconds in F555W and 1,200 seconds 
in F814W. Our exposure times are about 40% larger in V and 50% larger in /. It is why, with the use of 
different methods/softwares in the photometric analysis, we can reach stars 0.5 mag and 1 mag fainter in V 
and /, respectively. 

These two CMDs show a relatively shallow RGB and a horizontal branch (HB) populated predominantly 
redward of the RR Lyrae instability strip. Both of these features suggest that Gl is a rather metal-rich 
stellar system. However, the CMD also reveals a blueward extension to the red HB clump composed of a 
small number of stars. All three of these populations were also noted by Rich etal. (1996) in their CMD of 
Gl. In particular, Rich etal. (1996) pointed out that the blue extension to the HB could possibly be the 
result of a chemical abundance spread in Gl. Indeed, the RGB does display a potentially significant color 
width. The statistical significance of this width is addressed below. 

We note that the morphology of Gl HB is more reminiscent of the HB morphology observed in the 
dwarf spheroidal galaxy Andromeda I (Da Costa etal. 1996) than in the globular cluster 47Tucanae 
(Vazdekis etal. 2001). This fits some of the conclusions of this work, which unveiled some characteristics of 
the stellar population of Gl, more typical of dwarf galaxies than of our classical view of globular clusters. 

In order to measure the magnitude of the HB, we construct a luminosity function of the data and 
fit a Gaussian curve to the most prominent peak in this luminosity function. This procedure yields 
V{HB) = 25.34 ± 0.07. The quoted error is the result of adding, in quadrature, estimated errors of ± 0.05 
in the determination of V(HB) and ± 0.05 in the photometric zeropoint. This is in excellent agreement 
with the work of Rich etal. (1996) who obtained V(HB) = 25.32 ± 0.05. To estimate the metallicity of 
Gl in a way that is independent of the photometric zeropoint and the reddening, we rely upon the slope of 
the RGB as calibrated by Sarajedini etal. (2000). Utilizing their measurement technique and calibration 
leads to a value of [Fc/H] = —0.95 ± 0.09 for the mean metal abundance of Gl on the scale of Zinn & 
West (1984). This abundance lies between the results of Rich etal. (1996) who obtained [Fe/H] ~ -0.7 on 
the scale of Zinn & West (1984), a value close to that of 47Tucanae, and those of Bonoli etal. (1987) and 
Brodie & Huchra (1990) who obtained [Fe/H] ~ -1.2. In addition, our abundance is in accord with the 
estimate of Stephens ct al. (2001) based on the V - K(RGB) of Gl; they find [Fe/H] = -0.9 ± 0.2. Lastly, 
we can utilize the mean RGB color along with the calculated metal abundance and Equ. 1 of Sarajedini 
(1994) to compute the reddening of Gl; we find E(V - I) = 0.05 ± 0.02. 

Let us return now to the RGB color width apparent in the CMD of Gl (see Fig. 2). If this feature 
is significant, i.e. not caused entirely by the photometric errors, then we can argue strongly that there is 
a metallicity dispersion in Gl. One method used to test this is to conduct artificial star experiments in 
order to estimate the true measurement error. To begin with, we select stars located along an RGB fiducial 
sequence. For each of two trials, we select 210 stars along this fiducial sequence and place them on the 
original PCI frames under the constraint that no two artificial stars be within two PSF radii of each other. 
The resultant images are then reduced with the same procedure as the original PC images. 

The filled circles in Fig. 3 represent the original magnitudes and colors of the 420 artificial stars while 
the open circles are their recovered values. We are interested in the color width of the artificial stars and 
how this compares with the observed RGB color width. The open circles in the lower panel of Fig. 3 show 
the color histogram of stars located 1.8 ± 0.25 magnitudes above the HB of Gl. This location is chosen 
because it minimizes the influence of asymptotic giant branch stars (Da Costa & Armandroff 1990; Geisler 
& Sarajedini 1999). The filled circles represent the color histogram of the artificial stars around the fiducial 
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sequence. Gaussian fits to these distributions yield a b s = 0.144 ± 0.010 for the observed width and 
(j err = 0.037 ± 0.004 for the width due to photometric errors. Subtracting these quantities in quadrature 
gives ai n t — 0.139 ± 0.011 for the intrinsic width of the RGB. We note in passing that these artificial 
star experiments do not provide a complete and total assessment of the photometric errors. Other sources 
of error, e.g., such as residual flat-field non- uniformities and residual dark current, are not included in our 
photometric error estimates. 

If we assume then that the intrinsic photometric width we calculate is due entirely to a metallicity 
dispersion in Gl, what is the corresponding range in [Fe/H] ? To estimate this quantity, we turn to the 
standard RGB sequences of Da Costa & Armandroff (1990), which we use to construct a relation between 
[Fe/H] and (V — /)o at 1.8 mags above the HB of the six standard clusters in that study. The resultant 
relation is quadratic, which means that the [Fe/H] range it implies for Gl depends on the reddening we 
assume. The relation is: 

[Fe/H] = -18.648 + 23.846(V - I) - 7.906(V - 1)1 (1) 

with (V - I) = 1.39 at V(HB) - 1.8. For example, if E(V - I) = 0.10, then we infer a 1-er [Fe/H] 
dispersion of crpc/H] = ± 0.50 dex; whereas if E(V — I) = 0.05, then crpe/H] = ± 0.39 dex. In any case, 
the metallicity dispersion in Gl is genuine and significant. In contrast, we applied the above technique to 
HST/WFPC2 photometry of the M31 globular cluster G219 (Ncill 2001). These images were reduced in the 
same manner as those of Gl presented herein. We find that the color width of its RGB is fully consistent 
with the photometric errors as expected for a system with a negligible metallicity dispersion. 

The fact that in their HST/NICMOS study of Gl, Stephens ctal. (2001) do not observe any spread in 
metallicity does neither contradict nor infirm our present result. This for two reasons. First, for any spread 
in [Fc/H], the corresponding spread in color is twice as small in the infrared (J — K) than in the visible 
(V — I). Second, their very small sample (they have about 200 stars while we have about 5,000 stars) would 
certainly affect statistically their measurement of any spread in metallicity. 

The only other globular cluster known to exhibit a significant metal abundance range is lu Centauri, 
the giant Galactic globular cluster (see Jurcsik 1998 for a compilation of abundance measurements). Its 
range in [Fe/H] is roughly one dex (Norris & Da Costa 1995), which is quite similar to the range inferred 
by our estimates for Gl. 

5. Ellipticity, Position Angle, and Surface Brightness Profile 

Surface photometry of the cluster is obtained from all available WFPC2 images, using the techniques 
described by Djorgovski (1988). We then combine the profiles extracted from different images, using the 
shorter, unsaturated exposures for the central part of the profile, and longer, higher S/N exposures (in 
which the cluster core was saturated) for the outer parts of the profiles. 

Table 2 gives the ellipticity e = 1 — b/a and the position angle PA as a function of the semi-major axis 
a, using the stack of all our F555W (V) frames along with short exposures obtained in the same filter and 
available in the STScI/HST archives. These data are displayed in Fig. 4. The ellipticity varies significantly 
with the semi- major axis a, from values smaller than e = 0.1 in the innermost and outermost parts of 
the cluster, but with values larger than e = 0.2 between 0.7" and 7", reaching a maximum of e = 0.3 at 
a = 2.1". The mean ellipticity of Gl is e ~ 0.2. The position angle PA is not significantly variable for 
semi-major axis values a larger than 0.2". There is no significant evidence for twist of isophotes. 
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Table 3 gives the surface brightness /Lty and integrated V magnitude as a function of the radius. The 72 
points of this observed surface brightness profile are displayed in Fig. 5. There is no observational evidence 
of the presence of unbound stars and/or tidal tails surrounding Gl, simply by the mere fact that we would 
need to reach stars a few magnitudes fainter than the turnoff to have a statistically significant sample of 
such escaping stars. 

In order to convert the observed count rates to the V band magnitudes, we used the standard 
transformations from Holtzman ct al. (1995), assuming for the color of the cluster (B — V) = 0.84 mag. 
Integration of the combined profile gives the total magnitude for the cluster, V int = 13.48 ± 0.05 mag, where 
the net estimated zero-point uncertainty includes both the errors due to the background subtraction, and 
the color transformation (they are of a comparable magnitude). This is in an excellent agreement with the 
values from van den Bergh (1969) and Reed et al. (1994), after applying the appropriate aperture and color 
corrections, which give 13.54 and 13.58 mag, respectively. We note that these ground-based measurements 
are likely to underestimate slightly the cluster brightness, due to the removal of the light covered by the 
bright foreground stars, which accounts for some of the systematic difference here. (None of these numbers 
include the Galactic extinction corrections.) 

6. Mass Estimators 

We have at our disposal two essential observational constraints allowing the mass determination of Gl. 
First, our new determination of its surface brightness profile (see Table 3 and Fig. 5) from HST/WFPC2 
images, providing essential structural parameters: the core radius r c = 0.14" = 0.52 pc, the half-mass 
radius = 3.7" = 14 pc, the tidal radius r t ~ 54" = 200 pc, implying a concentration c = log (r t /r c ) ~ 
2.5. Second, our already published central velocity dispersion from KECK/HIRES spectra, providing an 
observed velocity dispersion a b s = 25.1 kms -1 , and an aperture-corrected central velocity dispersion a p (0) 
= 27.8 kms" 1 (Djorgovski etal. 1997). 

6.1. King model and Virial mass estimates 

Masses of dynamical systems are difficult to evaluate, with different methods providing rather different 
results, and the scatter between the different values is generally much larger than their formal uncertainties. 
Consequently, it is worth presenting results from more than one method, thus giving a realistic estimate 
of the true uncertainty. We can first obtain simple mass estimates from King models and from the Virial 
theorem (see, e.g., Illingworth 1976). 

The first estimate, the King mass M-k, is given by the simple equation: 

M K = Pcr\[i = 167 r cl iu 2 (2) 

where the core radius r c = 0.52 pc, the dimensionless quantity /i = 220 for c = log (r t /r c ) = 2.5 (King 
1966), and the central velocity dispersion ct p (0) = 27.8 kms -1 . These values determine a total mass for the 
cluster of Mtot = 15 x 1O 6 A^0 with the corresponding M/Ly — 7.5 in solar units. 

The second estimate, the Virial mass Myir, is given by the simple equation: 



Mvir = 670 r h a 2 



(3) 
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where the half-mass radius = 14 pc and central velocity dispersion cr p (0) = 27.8 kms -1 . These values 
determine a total mass for the cluster of Aitot — 7.3 x 10 6 .M© with the corresponding M.jLy ~ 3.6 in 
solar units. The internal error of each of these mass estimates amounts to about 10%. 

The large difference between these two estimates is not particular to the present cluster, but due to the 
idiosyncrasies of each method and typical of these two mass estimators applied to any globular cluster. See, 
e.g., Table 6 showing the results of similar mass estimates in the case of u> Centauri, the brightest and most 
massive Galactic globular cluster. 

7. King-Michie model mass estimates 

The two above observational constraints, viz., surface brightness profile and central velocity dispersion, 
allow the use of a multi-mass King-Michie model as defined by Gunn & Griffin (1979). See Meylan (1987) 
and Meylan etal. (1994, 1995) for the case of u Centauri. 

7.1. The model 

The "lowered Maxwellian" energy dependence [exp(— A^E) — 1] has been shown (King 1966) to be a 
good approximation to the solution of the Fokker-Planck equation describing the phase-space diffusion and 
evaporation of stellar systems like globular clusters. Following Eddington (1916), Michie (1963) introduced 
possible radial (u 2 ^ Vg = anisotropy by a factor of the form exp(— /3J 2 ), where J is the total angular 
momentum. Da Costa & Freeman (1976) showed that single-mass, isotropic King models are unable to 
fit the entire surface brightness profile of most globular clusters; they generalized these simple models to 
produce more realistic multi-mass models with full equipartition of energy in the centre. 

In the present work, we use a multimass anisotropic King-Michie dynamical model, based on an 
assumed form for the phase-space distribution function in an approach similar to that of Gunn & Griffin 
(1979). Each of the twelve subpopulations used has an energy-angular momentum distribution function 
fi(E, J) such that: 

fi(E, J) cx [cM-^E) - 1] exp(-/3J 2 ). (4) 

Thermal equilibrium is assumed in the cluster center, because of the short expected relaxation time, in 
order to force Ai to be proportional to the mean mass m» of the stars in the subpopulation considered. A 
model is specified by an initial mass function (IMF) exponent x (see below) and by four parameters: (i) the 
core radius r c , (ii) the scale velocity v s , (iii) the central value of the gravitational potential W , and (iv) the 
anisotropy radius r a , beyond which the velocity dispersion tensor becomes increasingly radial. 

7.2. The initial mass function 

There is no observational determination of the present-day mass function in Gl. Consequently, and in 
order to mimic a real cluster, main sequence (MS) stars, white dwarfs (wd) and heavier remnants (hr), such 
as neutron stars, are distributed into twelve different mass classes, adopting the usual power-law form for 
the initial mass function: 

dN cx m^ x d\og(m) (5) 



- 9 - 



where the exponent x would equal 1.35 in the case of Salpeter (1955). 

This initial mass function must be cut off at both extremities. The upper limit has neither dynamical 
nor photometric influence because of the rather small initial number of massive stars which have in any 
case already evolved into dark stellar remnants. This upper cutoff is chosen rather arbitrarily at 100 Mq. 
The lower limit is much more controversial because of the potential dynamical importance of numerous 
low-mass stars with low-luminosity. As found by Gunn & Griffin (1979), this lower mass cutoff, if it is low 
enough, does not significantly affect the cluster structure as traced by the giant stars. Numerous light stars 
do not change the quality of the fit, they simply increase the cluster mass and mass-luminosity ratio. The 
individual mass of the lightest stars is taken equal to 0.13 Mq. 

Owing to the total lack of observational constraints on the present-day mass function along the main 
sequence, the exponent x of the initial mass function is allowed to take different values in the following 
three mass intervals: Xhr, describing the heavy stellar remnants, resulting from the already evolved stars 
with initial masses in the mass range between 0.88 and 100 M Q ; x^ s , describing the stars still on the main 
sequence, with initial masses in the mass range between 0.32 and 0.88 M Q ; and x l ™ s describing the stars 
still on the main sequence, with initial masses in the mass range between 0.13 and 0.32M Q . 



7.3. The fit 

Both the HST/WFPC2 images providing the surface brightness profile (Table 3 and Fig. 5), and the 
integrated light KECK/HIRES spectra providing the central velocity dispersion (Djorgovski ctal. 1997) are 
heavily dominated by the light emitted by the brightest stars in Gl. All of these stars, viz., the giants 
and subgiants visible in the CMD (see Fig. 2), have individual masses equal to or slightly smaller than the 
turnoff mass, and belong to the same subpopulation containing stars with individual masses in the range 
0.63 to 0.88 Mq. Consequently, the fits between the models and the observations are made by using only 
this subpopulation, which contains the brightest members in the cluster, i.e., the giants, subgiants, turnoff 
stars, and the stars at the top of the main sequence. An acceptable model must match, simultaneously, 
first, the observed surface brightness profile, fitted by least squares to the integrated density profile of 
the subpopulation containing the bright stars, as determined by the model (see Fig. 5), and, second, the 
observed value of the velocity dispersion in the core of the cluster, which is unfortunately less constraining 
than the full velocity dispersion profile, as available, e.g., in the case of wCentauri (Meylan etal. 1995). In 
addition to the above two requirements, a model has to recover the total luminosity of the cluster, viz., My 
= -10.86 mag, to within 0.1 mag in order to be considered as satisfying. 



7.4. Relaxation time 

Two different "average" relaxation times are obtained for each model: a half-mass relaxation time and 
a central relaxation time. The term "average" means that these times depend on the mean stellar mass of 
the system, instead of being related to one particular species. Spitzer & Hart's (1971) standard formula 
transforms into: 

t rh = (8-92 10V) x x - {r * /lP K f: - (6) 

K y ' (M/lm Q ) log(0 AM/m) V ' 



where m is the mean stellar mass of all the stars in the cluster, M. the total mass of the cluster, and rt the 
half-mass radius. From Lightman & Shapiro (1978), the central relaxation time is given by: 



t r , = (1.55 10V) x 



(lm /m o ) v s 




(7) 



log(0.5M/m) lkm/ s 



where m is the mean mass of all the particles in thermal equilibrium in the central parts, v s the velocity 
scale, and r c the core radius. 



An extensive grid of about 150,000 models is computed in order to explore the parameter space defined 
by the Initial Mass Function (IMF) exponents x (defined over three different mass ranges Xhr, and 
x l Mg, and where x would equal 1.35 in the case of Salpeter 1955), the central gravitational potential W , 
and the anisotropy radius r a . 

Good models are considered as such not only on the basis of the x 2 of the surface brightness fit - 
the topology of the x 2 surface has no unique minimum — but also on the basis of their predictions of the 
integrated luminosity and mass-to-light of the cluster. We present hereafter results for the twelve models 
with lowest chi-square and fulfilling the other two requirements. 

The different columns in Table 4 give, for each model defined by its index, the central value of the 
gravitational potential W ; its MS mass function exponents x u ^ s and x l ^" s ; the fraction of its total 

mass in the form of stellar remnants such as neutron stars and white dwarfs; its concentration c = log 
{ft/r c ), core radius r c , half- mass radius r^, and tidal radius r t ; its central value of the mass density p Cl 
mean mass density ph inside the half-mass radius, and mean mass density pt inside the tidal radius. 

The different columns in Tabic 5 give, for each model defined by its index, the total mass Mtot of the 
cluster, in millions of solar masses, and its corresponding mass-to- light ratio M/Ly, also in solar units; the 
half-mass relaxation time t r h from Eq. (6), and central relaxation time t rj0 from Eq. (7). 

Since the velocity dispersion profile is reduced to one single value — the central velocity dispersion — 
the models are not constrained as strongly as in the case of u> Centauri (Meylan ctal. 1995), and equally 
good fits are obtained for rather different sets of parameters. Consequently, reliable results only relate 
to general parameters like the concentration and the total mass, but probably fail in any more detailed 
parameters. Nevertheless, very large areas of the parameter space can be eliminated with confidence since 
they never generate any successful models. 

The IMF exponent Xhr, describing the amount of heavy stellar remnants, appears in all good models 
to be very close to x — 1.35 (Salpeter 1955). The IMF exponent x^ s , describing the upper part of the MS, 
appears in all good models to be very close to x = 1.55. The IMF exponent x l £fg, describing the lower part 
of the MS, is not very well constrained; this is in agreement with Gunn & Griffin's (1979) findings that the 
lower-mass stars do not significantly affect the cluster structure as traced by the giant stars. The fraction 
of the total mass of the cluster in the form of heavy stellar remnants (neutron stars and white dwarfs) is 
always in the range of 18 to 20%. 

With a concentration c = log {r t /r c ) somewhere around 2.5, Gl is definitely a very concentrated 
globular cluster: all good models present clearly all the characteristics of a collapsed cluster. With its small 
core radius, Gl is hardly resolved in its central parts while its envelope exhibits an extended profile typical 



7.5. 



Results from King-Michie models 
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of a collapsed cluster (see also Djorgovski 1988). The core radius has a mean value of about r c ~ 0.52 pc, 
the half-mass radius ~ 13.5 pc, and the tidal radius, the most uncertain of these three radii, has a mean 
value of about r t ~ 200 pc. The corresponding mass densities have mean values of about p c = 4.7 x 10 5 
M Q pc- 3 , p h = 7.5 x 10 2 M pc- 3 , p t = 4.2 x lO" 1 M Q pc~ 3 . 

With a total mass somewhere between 14 and 17 10 6 Mq, and with the corresponding mass-to- light 
ratio A4/Ly between 6.9 and 8.1, Gl is significantly more massive than u) Centauri, maybe by up to a factor 
of three. These King-Michie mass estimates are in full agreement with the King mass estimate (A4k = 15 
x 1O 6 X with M/L v ~ 7.5), while the Virial mass estimate {Mvir = 7.3 x 10 6 A4© with M/L v ~ 3.6) 
is smaller by about a factor of two. It is worth mentioning that such a difference between King and Virial 
mass estimates is not specific to Gl: the same factor of about two is also observed between the King-Michie 
and Virial mass estimates of any cluster. See, e.g., in Table 6 the comparison of the results obtained for Gl 
and u Centauri (Meylan & Mayor 1986, Meylan 1987, Meylan etal. 1995, and this paper). 

8. Is Mayall II = Gl a genuine globular cluster ? 

We reach the following conclusions: (i) Gl is only the second globular cluster in which we observe 
a significant metallicity spread among its giant stars, u> Centauri being the first such case; (ii) All mass 
estimates (Table 6) give a total mass for Gl equal to as much as three times the total mass of u) Centauri; 
(iii) With c = log (r t /r c ) ~ 2.5, Gl is significantly more concentrated than 47Tucanae, which is a massive 
Galactic globular cluster considered on the verge of core collapse; all Gl structural parameters deduced 
from its surface brightness profile as well as the model densities are typical of a collapsed cluster; (iv) Gl 
is the heaviest of the weighed globular clusters. 

Although lo Centauri is by far the brightest and most massive globular cluster in our Galaxy, Gl may 
not be the only such massive globular cluster belonging to M31. This galaxy, which has about three times 
as many globular clusters as our Galaxy, has at least three other bright clusters which have central velocity 
dispersions larger than 20 kms -1 (Djorgovski etal. 1997). Unfortunately, so far, Gl is the only such cluster 
imaged with the high spatial resolution of the HST/WFPC2 camera, and consequently the only such 
massive cluster with a high quality surface brightness profile and known structural parameters. Gl and the 
other three bright M31 globular clusters represent probably the high-mass and high- luminosity tails of the 
mass and luminosity distributions of the rich population of M31 globular clusters. 

The above large mass estimates, implying a rather deep potential well, obviously relate to the 
metallicity spread whose origin is still unknown. There are essentially three different scenarios to explain 
such a metallicity spread: (i) metallicity self-enrichment in the globular cluster, (ii) primordial metallicity 
inhomogencity in the proto-cluster cloud, and (iii) the present globular cluster is merely the remaining core 
of a previously larger entity. 

Even more so than u) Centauri, Gl could be a kind of transition step between globular clusters and 
dwarf elliptical galaxies, in being the remaining core of a dwarf galaxy whose envelope would have been 
severely pruned by tidal shocking due to the bulge and disk of its host galaxy, M31. 

Kormendy (1985, 1987, 1994) used the four following quantities — the central surface brightness 
p,(0,V), the central velocity dispersion er p (0), the core radius r c , and the total absolute magnitude My - 
in order to define various planes from combinations of two of the above four parameters, e.g., (p(0,V) vs. 
log r c ). In all four planes plotted by Kormendy (1985, his Fig. 3), the various stellar systems segregate into 
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three well separated sequences: (i) ellipticals and bulges, (ii) dwarf ellipticals, and (iii) globular clusters. 
When plotted on any of the four planes, Gl appears always on the sequence of globular clusters, and 
cannot be confused or assimilated with cither ellipticals and bulges or dwarf ellipticals. The same is true 
for lo Centauri. Consequently, in spite of their large masses and internal stellar metallicity spreads, Gl and 
oj Centauri look like genuine bright and massive globular clusters. But where on these diagrams would the 
remaining cores of dwarf galaxies be located ? 

Because of its very large central velocity dispersion, M32 could not be the progenitor of a star cluster 
like Gl (van der Marel etal. 1998). But this is not true for the nucleated dwarf galaxy NGC 205, which has 
a central velocity dispersion similar to Gl (Carter & Sadler 1990, Held etal. 1990, Bender etal. 1991). The 
nucleus of NGC 205 is the only one for which the values of the four parameters used in Kormendy's diagrams 
are known, viz., the central velocity dispersion ct p (0) = 30 kms -1 (Bender etal. 1991), the central surface 
brightness /z(0,V) = 12.84 mag arcsec -2 , the core radius r c — 0.35 pc, and the total absolute magnitude 
My = -9.6 (Jones etal. 1996). These values place NGC 205 nucleus, in the Kormendy's diagrams, very 
close to Gl, right on the sequence of the globular clusters. Although this does not prove that all (massive) 
globular clusters are the remnant cores of nucleated dwarf galaxies, it shows that at least the nucleus of one 
such dwarf exhibits characteristics identical to those of globular clusters. It would be useful to know more 
about the nuclei of other dwarf galaxies. 

Of these massive globular clusters, the most nearby, ui Centauri, is naturally the best studied, 
nevertheless without decreasing the number of its conundrums. For instance, recent photometric (Anderson 
1997) and kinematic (Norris et al. 1997) studies show that u Centauri presents numerous characteristics 
about its stellar populations which remain without any explanation, if not completely puzzling. Presently, 
only the anti-correlation of metallicity with age (Hughes & Wallerstein 2000) and the unusual patterns of 
CN elements (Hilker & Richtler 2000) recently observed in ui Centauri suggest] that this cluster enriched 
itself over a timescale of about 3 Gyr. This contradicts the general idea that all the stars in a globular 
cluster are coeval, and may favor the origin of w Centauri as being the remaining core of a larger entity, e.g., 
of a former nucleated dwarf elliptical galaxy. Such a general idea had already been suggested by Zinnecker 
(1988) and Freeman (1993). In any case, the very massive globular clusters, by the mere fact that their 
large masses imply complicated stellar and dynamical evolution, may blur the former clear (or simplistic) 
difference between globular clusters and dwarf galaxies. 
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94A, GO-05907.01-94A, and GO-06477.02-95A from the Space Telescope Science Institute, which is operated 
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SGD was supported in part by the grant GO-060399 from STScI. GM thanks Jennifer Lotz (JHU), Roeland 
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Table 1. General information about Mayall II = Gl 



Parameters Mayall II = Gl 

a Gl (J2000) 00 h 32 m 46.6 s 

<5 Gl (J2000) +39° 34' 40" 

a M31 (J2000) 00 h 42 m 44.5 s 

<5 M31 (J2000) +41° 16' 29" 

Distance D to M31 770 kpc 

Color excess E(B — V) 0.06 mag 

True distance modulus (m — M) 24.42 mag 

Observed magnitude V 13.48 mag 

Absolute magnitude My —10.94 mag 

Central V surf bright fx(0,V) 13.47 mag arcsec~ 2 

Age ~ 15 Gyr 

Metallicity [Fc/H] -0.95 

Mean cllipticity e 0.2 

Radial velocity V r -331 ± 24 km s" 1 

Velocity dispersion (T (, s 25.1 kms -1 

Vel. dis. aperture corrected <r(0) 27.8 kms -1 
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Table 2. Mayall II = Gl: ellipticity e and position angle PA as a function of the semi-major axis a 
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Table 3. Mayall II = Gl: surface brightness [iy and integrated V magnitude as a function of the radius 
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14.767 


10.002 


22 


.768 


± 





.226 


13.543 


0.598 


15 


.808 


± 


0.050 


14.752 


10.226 


22 


.784 


± 





.225 


13.540 


0.638 


15 


.926 


± 


0.051 


14.693 


11.179 


23 


.063 


± 





.266 


13.528 


0.710 


16 


.133 


± 


0.055 


14.600 


12.084 


23 


.283 


± 





.320 


13.518 


0.763 


16 


,287 


± 


0.058 


14.541 


13.919 


23 


.849 


± 





.458 


13.502 


0.858 


16 


.580 


± 


0.066 


14.452 


14.269 


23 


.991 


± 





.516 


13.500 


0.868 


16 


.577 


± 


0.064 


14.444 


14.598 


24 


.113 


± 





.556 


13.498 


0.974 


16 


.836 


± 


0.065 


14.363 


17.632 


2d 


.941 


± 





.987 


13.484 


1.036 


16 


.940 


± 


0.064 


14.322 


18.214 


25 


.051 


± 


1 


.054 


13.482 


1.182 


17 


.188 


± 


0.057 


14.237 


18.945 


25 


.190 


± 


1 


.137 


13.480 
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Tabic 4. Mayal II = Gl: multi-mass King-Michie model parameters 



model 


Wo 


up 
X MS 


low 
X MS 




c 


r c 


rh 


n 


po 


Ph 


Pi 


index 








% 


log (rt/r c ) 


[PC] 


[PC] 




[M oP c- 3 ] 


[M pc- 3 ] 


[M pc- 3 ] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


1 


15.25 


1.50 


-0.50 


20 


2.49 


0.53 


12.3 


163 


4.7E+5 


8.9E+2 


7.8E-1 


2 


15.50 


1.55 


-0.50 


20 


2.59 


0.52 


13.2 


199 


4.9E+5 


7.7E+2 


4.5E-1 


3 


15.50 


1.55 


-0.40 


19 


2.57 


0.52 


13.2 


193 


4.8E+5 


7.7E+2 


4.9E-1 


4 


15.50 


1.55 


-0.30 


19 


2.55 


0.53 


13.2 


187 


4.7E+5 


7.8E+2 


5.4E-1 


5 


15.75 


1.55 


0.00 


18 


2.61 


0.52 


13.9 


212 


4.8E+5 


6.9E+2 


3.9E-1 


6 


15.50 


1.60 


-0.50 


19 


2.56 


0.53 


13.2 


193 


4.7E+5 


7.7E+2 


4.9E-1 


7 


15.50 


1.60 


-0.40 


19 


2.55 


0.53 


13.2 


187 


4.6E+5 


7.7E+2 


5.4E-1 


8 


15.75 


1.60 


-0.30 


19 


2.65 


0.51 


14.1 


230 


4.9E+5 


6.7E+2 


3.1E-1 


9 


15.75 


1.60 


-0.20 


18 


2.63 


0.52 


14.0 


221 


4.8E+5 


6.8E+2 


3.5E-1 


10 


15.75 


1.60 


-0.10 


18 


2.61 


0.52 


14.0 


213 


4.7E+5 


6.8E+2 


3.9E-1 


11 


15.75 


1.60 


0.00 


18 


2.59 


0.53 


14.0 


205 


4.6E+5 


6.9E+2 


4.4E-1 


12 


16.00 


1.60 


+0.20 


17 


2.68 


0.52 


15.0 


245 


4.8E+5 


5.9E+2 


2.6E-1 



Table 5. Mayal II = Gl: multi-mass King-Michie model parameters 



model 


Mtot 


M/L v 




tj-,o 


index 


[1O 6 .M ] 


[©] 


[10 9 yr] 


[10 6 yr] 


(1) 


(2) 


(3) 


(4) 


(5) 


1 


14 


6.9 


43 


15 


2 


15 


7.1 


50 


14 


3 


15 


7.2 


50 


15 


4 


15 


7.3 


50 


15 


5 


16 


7.6 


57 


14 


6 


15 


7.3 


51 


15 


7 


15 


7.3 


51 


15 


8 


16 


7.5 


57 


14 


9 


16 


7.6 


57 


15 


10 


16 


7.7 


58 


15 


11 


16 


7.8 


58 


15 


12 


17 


8.1 


67 


15 



Table 6. Three different mass estimates for the two globular clusters Mayall II = Gl and NGC 5139 = 

uu Ccntauri 



Mass 


Gl 


Gl 


lo Cen 


u) Cen 


Estimator 


mass [1O 6 .M0] 


M/L v le] 


mass [10 6 A4q] 


M/L v [e] 


King 


15. 


7.5 


4.3 


3.5 


Virial 


7.3 


3.6 


3.0 


2.4 


King-Michie 


14-17 


6.9-8.1 


5.1 


4.1 
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Fig. 1. — Mayall II = Gl as seen with the PC chip of the WFPC2 camera onboard the Hubble Space 
Telescope. The cluster is surrounded by two bright foreground stars. This image is a composite from F555W 
(V) and F814W (I) frames with a total exposure time of 2,200 seconds in V and of 1,800 seconds in /. The 
held is 31.5" x 31.5". North is 131° clockwise from vertical direction and East is 41° clockwise from vertical 
direction. (A non-degraded color version of this figure is available from the first author) . 



- 21 - 



CM 

IN P — 1 1 — 



-1 1 1 1 r- 



CM 



—i — i — i — i — i — i — i — i — i — i — i — i — i — i — 



CD 
CM 



V-l 



V-l 



Fig. 2. — Color-magnitude diagram (CMD) of 

Mayall II = Gl. The left panel shows the V vs. V ~ I CMD, while the right panel displays the I vs. V — I 
CMD. Each panel contains the same 4903 stars. 
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Fig. 3. — Artificial star experiments. Upper panel: color-magnitude diagram where filled circles represent 
the original magnitudes and colors of the 420 artificial stars while the open circles are their recovered values. 
The dots represent the genuine Gl stars. Lower panel: the color widths of artificial and observed stars. The 
open circles show the color histogram of stars located 1.8 ± 0.25 magnitudes above the horizontal branch of 
Gl while filled circles show the color histogram of the artificial stars. 
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Fig. 4. — Ellipticity and position angle as a function of the semi-major axis, for Mayall II = Gl, using the 
stack of all F555W (V) frames obtained with the PC chip of the WFPC2 camera onboard the Hubble Space 
Telescope. 



i — -24-' 




1 2 

log R [Rc] 

Fig. 5. — Surface brightness profile of the globular cluster Mayall II = Gl, from HST/WFPC2 shallow 
and deep images in F555W ~ V filter; the continuous line represents a King-Michic model (first model in 
Table 4) fitted to the observed profile (Meylan etal. 1999). 



